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Small-angle Bhabha scattering data recorded at the Z resonance and large-angle 
Qh| Bhabha scattering data recorded at ^Js = 189 GeV by the L3 detector at LEP are 

used to measure the running of the effective fine-structure constant for spacelike 
momentum transfers. The results are 

a'\-2.l GeV 2 ) - a-\-6.25 GeV 2 ) = 0.78 ± 0.26 
c3 ' a-V-12.25 GeV 2 ) - cT 1 (-3434 GeV 2 ) = 3.80 ± 1.29, 



in agreement with theoretical predictions. 
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Introduction 



At zero momentum transfer, the QED JTj fine structure constant a(0) is very accurately known 
from the measurement of the anomalous magnetic moment of the electron and from solid-state 
physics measurements j2j: 

a -1 (0) = 137.035 999 76 (50). 

In QED, vacuum polarization corrections to processes involving the exchange of virtual photons 
result in a Q 2 dependence, or running, of the effective fine-structure constant, a(Q 2 ). This Q 2 
dependence is usually parametrised [3 as 

Whereas the leptonic contributions to Aa(Q 2 ) can be accurately calculated, due to non- 
perturbative QCD effects the vacuum polarization contributions from quark loops cannot be 
calculated exactly. Therefore, dispersion integral techniques |IHEj are used to estimate these 
contributions from the measured cross sections for the process e + e~ — > hadrons at low centre- 
of-mass energies, yielding at Q 2 = m|: a _1 (m|) = 128.886 ± 0.090. Similar evaluations have 
been made in Reference [7j and, under additional theoretical assumptions, in References [HUH]. 

Processes involving photon exchange at non-zero momentum transfer Q 2 yield amplitudes 
proportional to a(Q 2 ). At e + e _ colliders this implies that the measurement of the relative rates 
of processes involving different Q 2 values gives access to the running of a(Q 2 ) between those 
Q 2 values. This was first exploited by the TOPAZ Collaboration ^U|, which derived a value 
of the fine structure constant at a scale of Q 2 = 3338 GeV 2 from the ratio of the muon pair 
annihilation cross section to the two-photon-induced muon pair cross section. 

In regions where Bhabha scattering is dominated by t-channel photon exchange, this process 
allows the study of a(Q 2 ) in the spacelike region, Q 2 = t = —s(l — cos6>)/2 < 0, where 9 is the 
angle of the outgoing e~ with respect to the e~ beam direction. The VENUS Collaboration 
has recently interpreted large-angle Bhabha scattering measurements in terms of the evolution 
of a(Q 2 ) in the range 100 GeV 2 < -Q 2 < 2916 GeV 2 P]. 

This article interprets measurements of Bhabha scattering by the L3 experiment at LEP in 
terms of the evolution of a(Q 2 ), using two complementary measurements: first, an analysis of 
the small-angle Bhabha scattering data in the polar angular range 32 mrad < 9 < 54 mrad, 
collected at centre-of-mass energies around the Z mass and used for the high-precision lumi- 
nosity measurement in the years 1993-1995; second, the measurement of the Bhabha scattering 
cross section in the polar angular range 20° < 9 < 36° performed at yfs = 188.7 GeV in 1998. 
They are used to study the running of a(Q 2 ) in the range 2.1 GeV 2 < — Q 2 < 6.25 GeV 2 and 
12.25 GeV 2 < -Q 2 < 3434 GeV 2 , respectively. 



Small-angle Bhabha Scattering 
Data Analysis 

The analysis of the small-angle Bhabha scattering data follows closely that of the luminosity 
measurement and makes use of the same detectors: 

• two small-angle calorimeters consisting of BGO crystals, providing a precise energy mea- 
surement for electromagnetic showers. The highest-energy cluster reconstructed in each 
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calorimeter is retained for analysis. The energy of one cluster must exceed 0.8£^be am ; the 
energy of the cluster in the other calorimeter must be greater than 0.4i?beam- To avoid 
edge effects their reconstructed polar angles must be well contained in the calorimeter, 
28 mrad < 9 < 65 mrad; 

• a silicon strip detector, consisting of two layers of r-measuring strips and one layer of 
^-measuring strips, installed in front of each BGO calorimeter. This detector allows a 
precise definition of the fiducial volume. The coordinates of the cluster reconstructed in 
the calorimeter are projected onto each r-measuring silicon detector plane in turn, and a 
matching window (5 mm for r-measuring strips and 2.5° for 0-measuring strips) is used 
to search for corresponding hit strips. If found, their coordinates are used. Otherwise the 
BGO coordinates are retained. 

Angular cuts are applied on the coordinate measurements for each r-measuring plane in turn. 
The fiducial volume is the same as that used for the luminosity measurement (32 mrad < 9 < 
54 mrad), divided into four polar angular bins, with boundaries at 32, 35, 40, 46, and 54 mrad. 
For each event, each of the four coordinate measurements (two on each side) is entered in the 
corresponding angular bin with a weight of 0.25. 

The silicon coordinate reconstruction is sensitive to malfunctioning silicon strips. Therefore, 
only those data are used for which the silicon detector was fully functional, resulting in a sample 
of 6.7- 10 6 Bhabha events corresponding to an integrated luminosity of 98.8 pb _1 . In addition, 
the polar angular bins are chosen to have boundaries sufficiently far away from malfunctioning 
strips (8 in 1993, 9 in 1994, and 11 in 1995) and from the edges of the flare in the beam pipe 
(42.5 and 50 mrad) on the outgoing e~ side. The data are grouped into 7 sets, according to 
the three years and the different centre-of-mass energies. 

Detector effects are corrected for, bin by bin, using a sample of 2.6 • 10 6 fully simulated ^Bj 
and reconstructed BHLUMI ^3] Monte Carlo events. The correction factors differ from 1 by 
at most 0.4% and have statistical uncertainties of about 4 • 10~ 4 . 

The data are compared with predictions from the BHLUMI Monte Carlo program. This 
program is modified to allow the running of a{Q 2 ) in the spacelike region to be different from 
the nominal one of Ref. [I]. The running in the timelike region is left unchanged. This is 
relevant for the interference term between the dominant t-channel photon-exchange and the 
s-channel Z-exchange diagrams, which contributes up to 0.15% to the Bhabha cross section in 
the luminosity monitor fiducial volume at centre-of-mass energies around the Z resonance. The 
predictions are calculated for each data set separately. 

Since the small-angle Bhabha scattering data are used for the absolute normalisation of all 
processes, only the shape of the polar angular distribution can be used in this study. Therefore, 
the compatibility of the data with any given parametrisation, p, of the running of a{Q 2 ) is 
evaluated in terms of a likelihood formula according to multinomial statistics: 



is the fractional cross section predicted in bin j (= 1 ... 4) for data set i {— 1 . . . 7) for parametri- 
sation p, rijj is the observed event weight in each angular bin, and P is the multinomial prob- 
ability distribution. 



ln£ = - ^lnP(n M ,...,n 4 ,;;/iV-->/T,i), 



(2) 



where 
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Results 



The data are first compared with two hypotheses: the theoretical prediction [I] (referred to 
as 'normal running') and the prediction for the case of no running of a{Q 2 ) between Q 2 and 
the lowest scale accessible in the analysis, Ql = —2.1 GeV 2 . Figure [T] shows, for the 1994 
data which have the highest statistics, the ratio of the measured and expected fractional 
event weights for the two hypotheses. The likelihood difference between the two hypothe- 
ses is — (ln£(no running) — ln£(normal running)) = 10.24 for the complete data set, excluding 
the hypothesis of no running by 4.5a. 

The data are then compared with a parametrisation in which a term, linear in Q 2 and with 
a slope, S, to be fitted, is added to the normal running. The Q 2 dependence of Equation [T] 
then becomes 

/ n 2N _ (o\ 

aW ] 1 - Aa(Q 2 ) - S -(Q 2 - Ql) ' {6) 

with Ql as given above. The fit yields the result S = (—3.6 ±2.7) • 10 -4 GeV -2 , where the error 
is statistical only. 

The following sources of systematic uncertainties on S have been studied: 

• Statistical uncertainties on the detector correction factors: their effect is studied by gen- 
erating a large number of random number sequences in which each correction factor is 
changed according to a Gaussian distribution having a width equal to its statistical un- 
certainty. Since the same Monte Carlo sample is always used for the corrections, the 
changes in the correction factors are assumed to be fully correlated between the same 
bins of different data sets. Each time the corrections are changed, and the spread in the 
fitted values of S is found to be 1.5 • 10 -4 GeV -2 ; 

• Statistical uncertainties on the theoretical predictions: their effect is studied in the same 
way. In this case, the statistical uncertainties on the cross sections are assumed to be 
fully correlated between different choices of the parametrisation of a(Q 2 ), but to be 
uncorrelated between different data sets. This yields an uncertainty of 0.5 ■ 10 -4 GeV -2 . 

• The size of the matching window used for the coordinate reconstruction: this is studied 
by changing it to 4 mm and 2° for the r- and 0-measuring strips, respectively, and the 
difference in the fit result (0.8 • 10 -4 GeV -2 ) is assigned as a systematic uncertainty. 

• The consistency between the results obtained using the individual r-measuring silicon 
layers: instead of combining the results from the four r-measuring layers, the analysis can 
also be performed for each side or layer separately. The fitted values of S are —6.7 ■ 10 -4 
using the layers on the outgoing e + (—z) side, and —0.5 • 10 -4 using the layers on the 
outgoing e - (+z) side. The results obtained using each of the two layers on either side 
are in excellent agreement with each other. 

From the ratio of the (corrected) fractions in each bin for the —z and +z sides, the 
difference is seen to be due to a rather poor overall agreement between the observed 
fractions for the two sides. The cause for this poor agreement is an incomplete simulation 
of the beampipe material traversed at small angles by the e 1 * 1 before entering the silicon 
detector. The uncertainty assigned to this source is half the difference between the results 
for the two sides, 3.1 • 10 -4 . This dominates the total measurement uncertainty. 

As a cross-check on the material effects, the last two bins are merged into a single bin 
which covers the complete region of increased material in the flare of the beam pipe. 
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Figure 1: Ratio of the measured and expected fractional event weights in each 
polar angular bin for the standard theoretical predictions (solid circles) and the 
assumption of no running (open circles), for the data collected in 1994. 
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An alternative is to take out from the analysis the first bin, which represents the most 
material for the layers on the other side. In both cases the results are stable within the 
systematic uncertainty assigned. 

Finally, the effect of the choice of polar angular bins is studied by performing the fits for 
several choices of binning. Variations of the size of the systematic uncertainty assigned 
are found, which are almost entirely due to variations observed on the +z side. Their 
cause is the same as that of the poor overall agreement, and no additional systematic 
uncertainty is assigned. 

• The uncertainty on the vertex position: this affects the two sides in opposite ways. Its 
effect has been estimated by artificially moving the vertex position, and is found to be 



The largest contribution (0.04%) to the theoretical systematic uncertainty on the luminosity 
measurement at centre-of-mass energies around the Z resonance [Tj)| comes from the theoretical 
uncertainty on the value of a(Q 2 ) in the small-angle Bhabha scattering region, and hence should 
not be considered in this analysis. Theoretical uncertainties on other contributions can affect 
the angular dependence only slightly. For example, the uncertainty due to missing 0(a s L 3 ) 
(where L represents the 'large logarithm' ln(|t|/m e )) contributions to the uncertainty on S is 
estimated to be about 4 • 10~ 6 , i.e. negligible compared to the experimental systematics. 
The final result of the measurement of S is thus 



Large-angle Bhabha Scattering 
Data Analysis 

The subdetectors used for the measurement of large-angle Bhabha scattering are the BGO 
electromagnetic calorimeter and the central tracker, described in detail in Reference JH]. The 
data were collected in 1998 at yfs = 188.7 GeV and correspond to an integrated luminosity of 
175.9 pb _1 . Electrons are identified as energy deposits of at least 0.5 GeV in the calorimeter, 
with an electromagnetic transverse shower shape, and at least six associated track hits within 
a three degree azimuthal angular range. At least one electron must be observed within the 
fiducial volume 20° < 9 < 36°, and one within 144° < 9 < 160°. At most six electromagnetic 
clusters are admitted in the electromagnetic calorimeter. The highest-energy electron must 
have an energy exceeding 0.5-Ebeam- In this angular range, t-channel photon exchange gives 
the dominant contribution to the cross section, with s — t interference and s-channel exchange 
contributing only -7% and 0.3%, respectively. 

The signal efficiency is estimated using 2.5 • 10 6 fully simulated and reconstructed BH- 
WIDE Monte Carlo events generated in the angular range 5° < 9 < 175°. The back- 
ground, consisting mainly of e + e~ — > t + t~(7) events, is small, 0.12%, and its uncertainty has 
a negligible effect on the cross section measurement. 

Given the large Q 2 range covered by this analysis, the parametrisation of Equation El is not 
appropriate. Deviations from the normal running of u(Q 2 ) are parametrised, for any Q 2 , using 
the following modification of Equation 21 



negligible. 



S = (-3.6 ± 2.7(stat.) ± 3.5(syst.)) • 10" 4 GeV" 2 . 



(4) 



a{Q 2 ) 



a(0) 



(5) 



l-C ■ Aa(Q 2 )' 
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where Aa(Q 2 ) is calculated according to Reference jlj. The evolution of a{Q 2 ) is then deter- 
mined from the ratio of measured to predicted cross sections, cr meas /a pTe( i(C), where a" P r e d(C) is 
the predicted cross section for a given value of C. Since Equation also modifies the value of 
a{Q 2 ) in the Q 2 range used for the luminosity measurement, also the measured luminosity and 
hence the measured cross section is changed for 1. This is accounted for by the use of the 
ratio 

n/y-A '-'meas ^lumi 

(c) 

H \ C ) = Z 77n ' ~ Tr 7v t b J 

Oprcdl^J ^lumi^ — 1 J 

where o" mmi (C) is the small-angle Bhabha scattering cross section in the fiducial volume used for 
the luminosity measurement. C is then determined by equating R(C) = 1. The cross section 
predictions are obtained using the BHWIDE and BHLUMI Monte Carlo programs for the large- 
angle and small-angle Bhabha scattering processes, respectively, modified to parametrise the 
evolution of a(Q 2 ) according to Equation 



Results 

The number of selected events is 23940. The selection efficiency is estimated to be (93.38 ± 
0.08)%, where the error is due to limited Monte Carlo statistics. The measured cross section is 
cr meas = 145.6 ± 0.9 pb. 

Systematic uncertainties on the cross section measurement are due to: 

• The description of the energy response: Figure 121 shows the energy of the highest-energy 
electron, normalised to -Ebeam- The low-energy tail of the distribution is not very well 
described by the Monte Carlo simulation. Its effect on the selection efficiency is estimated 
by adjusting the Monte Carlo distribution to obtain a better agreement, and is found to 
be negligible. 

• The definition of the fiducial volume: this is estimated by varying the angular cuts and 
results in a 0.4% uncertainty. 

• The cut on the number of track hits: a comparison between data and Monte Carlo events 
has been used to adjust the Monte-Carlo description of the single-hit efficiency of the 
inner tracker wires. Varying the cut on the number of hits changes the efficiency by 0.3%. 

Summing the individual contributions in quadrature, a total systematic uncertainty of 0.5% is 
obtained. The result of the cross section measurement is thus 

(T mea s = 145.6 ± 0.9(stat.) ± 0.8(syst.) pb. (7) 

The measured cross section agrees well with the theoretical prediction of a pre d(C = 1) = 
145.9 pb. Using Equation Inland assuming a theoretical uncertainty on o" pro d(C = 1) of 1.5% [T5] . 
the value of C obtained is 



C = 0.97 ± 0.12(stat.) ± 0.10(syst.) ± 0.29(theory). 



Interpretation of Results 

The results can be interpreted as a measurement of the evolution of a~ 1 (Q 2 ) in the spacelike 
region between the momentum transfer scales relevant for the analyses. In the small-angle 
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Figure 2: Distribution of the highest electron energy normalised to the beam energy. 
The arrow shows the position of the applied cut. 
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Bhabha scattering analysis these are taken to be the lowest and highest Q 2 values accessible, 
—2.1 GeV 2 and —6.25 GeV 2 . The difference is found to be 

a" 1 (-2.1 GeV 2 ) - a -1 (-6.25 GeV 2 ) = 0.78 ± 0.26, 

where the error reflects the total experimental uncertainty. 

In the large-angle Bhabha scattering analysis the relevant momentum transfer scales are 
taken to be the average Q 2 value used for the luminosity measurement, —12.25 GeV 2 , and 
the average Q 2 value for the t-channel contribution to the large-angle Bhabha scattering cross 
section, —3434 GeV 2 . The result is 

a~ 1 (-12.25 GeV 2 ) - oT 1 (-3434 GeV 2 ) = 3.80 ± 0.61(expt.) ± 1.14(theory). 

The results are displayed in Figure El For the purpose of the figure, for both measurements 
the value of a{Q 2 ) is fixed to its expectation at the lower momentum scale involved in the 
analysis. 

In conclusion, Bhabha scattering at LEP has been used to study the running of the fine- 
structure constant, a(Q 2 ), in the spacelike momentum transfer regions 2.1 GeV 2 < — Q 2 < 
6.25 GeV 2 and 12.25 GeV 2 < -Q 2 < 3434 GeV 2 . The data clearly establish a nonzero running 
as predicted by QED. 
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